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Abstract. Irreversibility is introduced to quantum graphs by coupling the graphs 
to a bath of harmonic oscillators. The interaction which is linear in the harmonic 
oscillator amplitudes is localized at the vertices. It is shown that for sufficiently strong 
coupling, the spectrum of the system admits a new continuum mode which exists 
even if the graph is compact, and a single harmonic oscillator is coupled to it. This 
mechanism is shown to imply that the quantum dynamics is irreversible. Moreover, it 
demonstrates the surprising result that irreversibility can be introduced by a "bath" 
which consists of a single harmonic oscillator. 



1. Introduction 

Quantum graphs emulate many properties of quantum chaotic systems, and their use 
is now widely spread in the physics and mathematics literature 121 El IH 03 EH El • In 
many potential applications the physical systems are irreversible, and it is desirable to 
include this feature to quantum graphs, and thus enlarge their range of applicability. 

The standard way to introduce irreversibility to a quantum system is by coupling 
it to a "bath" which consists of a set of "irrelevant" degrees of freedom. The bath 
dynamics is governed, however, by a well defined quantum hamiltonian [HI EH E3 ITT] . 
The enlarged system, composed of the original system and the bath, is a proper quantum 
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system, with hamiltonian: 

H _ H sy stern + jjbath + w ^ 

^system an( j j^bath are operators in the "system" and "bath" Hilbert spaces, respectively, 
and W is the coupling. In the physics literature, irreversibility is generally associated 
with the bath having a continuous spectrum. The physical justification is based on 
the observation that the Poincare recurrence time is infinite, or stated differently, the 
mean probability to return back to the initial state vanishes, which is a prerequisite for 
irreversibility. This argument can also be formulated mathematically. To simplify the 
presentation, assume that H system has a discrete spectrum, while the bath spectrum 
is continuous, with a spectral measure dfi(S) and 5 normalized eigenfunctions \xs)- 
Denoting by ty-p the projection of an eigenstate of H on the system subspace, The 
Schrodinger equation can be reduced to 

EV V = H s y stem V v + lim / dfi(£)W v (£) = I ■ . Wi{S)V v , (2) 

V\0 J Cj — C + 17] 

where Wp(£) is the projection of H^lx^:) on the system subspace. Using the identity 

lim — — = -mS{x) +V ( - ) , (3) 



r?\0 x + %r\ \x 
one traditionally associates irreversibility with the operator 

T(E) = -m J dn(£)6(E - £)W V (£)W},(£) . (4) 

T(E) is an anti hermitian operator, which introduces decay to the quantum time 
evolution. (This can be immediately shown when T(E) can be approximated by an 
energy independent operator). Two important points must be emphasized. First, the 
discussion presented above is valid only when the bath spectrum is continuous. Second, 
T(E) vanishes unless E is in the support of the spectral measure. Similar expressions 
appear e.g., for the width of resonances in the Breit-Wigner theory [T2*] . 

Various approximations and models where used to extract the main features of 
quantum irreversibility as described by (j3J). Amongst the most popular is the harmonic- 
bath model where the bath consists of a continuous set of harmonic oscillators which 
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are linearly coupled to the system degrees of freedom ^U]. In this case, some of the 
integrations can be carried out, and using first order perturbation theory, the decay 
width introduced by the bath can be computed explicitly 

The purpose of the present work is to carry out the above program for a system 
which consists of a graph that is coupled linearly to a harmonic bath. The system will be 
described in section (J2J) and the corresponding Schrodinger operator will be discussed. 
Because of the form of the coupling, the question of finding a self-adjoint extension is not 
simple. In the present work, a physically motivated form of the Schrodinger operator 
is constructed, and the rigorous proof of its validity is given in an adjacent paper by 
Michael Solomyak [T3*] . 

Irreversibility is demonstrated by computing the vertex scattering matrices Q3I2], 
which are responsible for the transfer of current between the bonds of the graph (section 
(jSJ). For sufficiently strong coupling, the vertex scattering matrices are sub-unitary, 
which means that the total probability current is not conserved. Rather, it decreases 
upon scattering at the vertex. In other words, the quantum mechanical evolution is 
not unitary. The surprise is that this happens even when a single harmonic oscillator 
is coupled at the vertex! This mechanism should not be confused with what happens 
naturally when an additional lead to infinity is coupled to the vertex, and its effect on 
the original scattering is to draw current. This mechanism which induces an "escape 
width" to the graph dynamics is not what is considered here. Rather, non unitarity is 
due to the coupling to a harmonic oscillator. 

The rest of the paper attempts to explain this result in physical terms. To this end, 
a simple graph is presented and solved approximately in section (J3J). Using this model, 
one can demonstrate and explain how the coupling of a single harmonic oscillator to 
a compact graph can modify the spectrum in a major way. If the coupling is week, 
the spectrum remains pure point. Beyond a critical strength, a continuous component 
appears in the spectrum. The appearance of the continuum occurs at the same value 
of the coupling constant at which the vertex scattering matrix becomes non unitary. A 



Irreversible quantum graphs 4 
rigorous treatment of this spectral problem is also provided in Solomyak's paper [T3] . 



2. The Schrodinger operator on a graph and the coupling to a harmonic 
bath 

A graph Q consists of V vertices connected by B bonds according to the connectivity 
matrix Cjj which takes the value 1 if the vertices i,j are connected, and it vanishes 
otherwise. We assign the natural metric to the bonds, and associate a length Lb to 
each. The position a; of a point on the graph is determined by specifying on which 
bond b it is, and its distance Xb from the vertex with the smaller index, < Xb < Lb- 
Sometimes, it will be convenient to denote the bond connecting the vertices % and j by 

(*,;)■ 

In the absence of a harmonic bath, the Schrodinger operator can be defined in the 
following way: Let x E Q and ^(x) a real valued and continuous function on Q, so that 
= ipb{%b) for x G b, and < Xb < Lb, where ipb{xb) are twice differentiable in 
the interior of the bond. We restrict ^(x) to be square integrable, and require that it 
is uniquely defined at the vertices. That is, all the functions ipb which correspond to 
bonds connected at a vertex i attain the same value <pi at the common vertex. Given 
an arbitrary set of V non-negative constants Aj one constructs the quadratic form: 

E A ^- ( 5 ) 



b=l 



(•Lb 

/ dx b 




h 


dxb 



i=l 



The Euler-Lagrange variational principle selects the stationary solutions of the quadratic 
form, as the solution of the Schrodinger equation on the bonds: 

~ E )Mx b ) = o, (6) 

subject to the boundary conditions 

- \4>i = o . (7) 



d 



The summation is over all the bonds (i, j) which emanate from the vertex i, and the 
derivatives are computed at the vertex. The quadratic form (0) is positive definite, and 
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therefore, the boundary condition (J7J) provides a self adjoint extension of the Schrodinger 
operator for any choice of the (non- negative) constants Aj. This operator determines 
the quantum dynamics of the "system". One can extend this operator to non compact 
graphs, by adding bonds from the vertices to infinity The corresponding scattering 
problem was amply discussed in the literature IH El EE Hj • 

The bath degrees of freedom consist of M harmonic oscillators with coordinates 
q m , and frequencies u m , subject to the hamiltonian 
1 / d 2 



Denote q = (<Zi, • • ■ , qu) and H osc = J2m=i ^m- H osc is the "bath" hamiltonian. The 
coupling of the bath to the system is introduced by extending the quadratic form (jSJ) to 
include the graph and harmonic bath coordinates. Consider *$?(x, q) = ipb{xb, q) for x G 
b, where each of the ipb{xb,<i) is square integrable. For every value of q, ^(x, q) is 
uniquely defined at the vertices, where, at the vertex i it assumes the value 0«(q). To 
each vertex, we assign a real valued function Aj(q), which is bounded at any finite 
domain of 1R M and is allowed to diverge algebraically as |q| — > oo. (In the following, we 
shall refer to these functions as the "form- factors" ). We construct the quadratic form, 

|d^ 6 (x 6 ,q) 2 



6=1 
V 



+ ipt(x b , q)H osc ^ b (x b , q )) 



dx b 

+ 5>(q)tf(q)} • (9) 

i=l 

Again, the Euler-Lagrange variational principle selects the stationary solutions of the 
quadratic form, as the solution of the Schrodinger equation on the bonds: 
d 2 \ 

^2 + H osc - EJ i[> b (x b , q) = , (10) 



subject to the boundary conditions 

- Ai(q)&(q) = . (11) 



v d 

^2 C ^aZ ^(m)( x (m)>^ 

j=1 ux (jj) 



The coupling between the "system" and the harmonic "bath" is thus mediated through 
the boundary conditions 
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When the form-factors A„(q) are non negative the quadratic form L osc [\l/] is positive 
definite, and the Schrodinger operator (fTUJl . subject to the boundary conditions (|TT|) is 
self adjoint. The standard linear coupling of harmonic baths is written as 

M 

Ai(q) = kim ' ( 12 ) 

m=l 

with A a,VxM matrix with q-independent real coefficients. Thus, L osc [\l/] is not definite, 
and the proof that the system (jlOlllj) is self-adjoint requires a special treatment, which 
is provided in The intuitive motivation behind this coupling scheme is that for 

a vertex with v = 2, and a single harmonic oscillator, the boundary condition at the 
vertex is equivalent to an effective potential Aq6(x). This represents a well localized 
coupling of the particle on the graph to the harmonic oscillator. 

3. The vertex scattering matrix 

The boundary conditions at the vertices of graphs (without coupling to any harmonic 
oscillators) (J7J) impose the conservation of the probability current across the vertex. In 
other words, the boundary conditions at a vertex can be described by a unitary vertex 
scattering matrix which provides the linear relation between the amplitudes of the waves 
which imping upon the vertex and the waves scattered from it p. The vertex scattering 
matrix which corresponds to the boundary conditions (J2J) can be written explicitly as 

-, . —2i arctan 
{ -\ A- e v i 

^j'(k) = + — (13) 

where i denotes the vertex under consideration, b and b' denote the bonds which emanate 
from the vertex, and k = \[E is the wave number under consideration. Note that this 
Vi x Vi vertex matrix distinguishes only between reflection (b = b') and transmission 
between different bonds (b ^ b'), which is independent of their specific identity. This is 
an expression of the fact that the boundary condition at the vertex is invariant under 
the interchange of the bonds. a^\,{k) is symmetric, and its unitarity can be easily 
demonstrated. The single vertex, and the bonds emanating from it, is some times 
referred to clS Qb "star-graph" . 
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Describing the quantum evolution and dynamics in terms of waves travelling freely 
along the bonds and scattered at the vertices, is a very useful concept in the theory 
of quantum graphs [T]. This is also true when harmonic oscillators are coupled, and 
we shall dedicate the rest of this section to the derivation of the the vertex scattering 
matrices which correspond to the boundary conditions We shall show, that the 

coupling of the harmonic oscillators can have a profound effect - the vertex scattering 
matrix is not necessarily unitary for large values of the coupling A. Thus, the evolution 
on the graph is not conserving and hence irreversible! 

To simplify the presentation, we consider the coupling to a single harmonic oscillator 
(M = 1), and choose the units such that u\ = 1. The spectrum of the combined system 
of a harmonic oscillator and the star graph without coupling, consists of overlapping 
continua with thresholds at E n = Further simplification is achieved by considering 

energies in the interval | < E < E\. In this energy range, and for x values away from 
the vertex, only the ground state (n = 0) of the harmonic oscillator is populated, which 
corresponds to purely elastic scattering. We shall compute the vertex scattering matrix 
by considering an incoming wave incident on the bond b' only. 



we expand the wave function in the complete orthonormal basis Xn(<?) of eigenf unctions 
of the harmonic oscillator hamiltonian h\{q). The bond wave functions are defined on 
the positive half line, with the common vertex at Xf, = 0. The expansion takes the form 



Denoting 



k 2 = E ; ko = k ; k n = \Jn — k 2 for n > 1 



(14) 




(15) 



The functions rjn (xb) are given by 




(16) 
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This way, ^(x, q) corresponds asymptotically to an incoming wave together with 
scattered waves in the ground state channel. In all other (closed) channels, the wave 
functions are evanescent along the bonds. Note that in (|16|) we made use of the expected 
invariance of the scattering matrix under the interchange of bonds, so that r stands for 
the reflection (diagonal) elements of the a matrix, and t stands for the transmission 
(off-diagonal) elements of a. Both parameters are independent of b. 

The continuity at of ty(x,q) at the vertex implies that the expansion coefficients 
dri ^ are the same for all 6, and therefore they will be denoted in what follows by a n . 
Moreover, 

l + r = t (17) 
Substituting in (fTTj) and using (fTTj) . we get 



2 A 1 

t = — + == ai, 

v 2% v/c /ci 



and, 



2«i + A( a 2 J -=- ) = ~^=t (19) 



k\ k 2 J \/k k 



Tb ~h 1 / Tl 

2a n + A ( a n+ i X l — Ya n _ x J-— ) = for n > 2 . 

Where, A = PJ) is an infinite set of linear, inhomogeneous equations, which are 

defined in terms of the infinite Jacobian matrix J(k, A) whose elements are explicitly 
given in (|19j). For large n 



l + 0[- \ . (20) 



We shall denote by Jo (A) the Jacobian matrix whose elements are defined by the 
equations with constant coefficients 

2a n + A (a n+ i + a n _i) = (21) 

which approximate 

Denoting by G{z\ k, A) the resolvent 
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we get 



Denote 



a\ = — lim G\ i (z'e; A:, A) — j t. (23) 

e/0 ' y/k fci 



-y(*,A) = ^r r limG lll (te;A ;) A) . (24) 

Substituting in IflB jl . we find 

2 (— 1 + -) + i7(/c, A) , 

1 — fy(k, A) 1 — Z7(fc, A) 

Note that when A = 0, the reflection and transmission coefficients reduce to (|13|) (with 
Aj = 0). The flux in the leads is conserved if (v — l)\t\ 2 + |r| 2 = 1, which is satisfied 
only if Im(7(fc, A)) = 0. We shall now show that this happens if and only if |A| < 1. In 
other words, the dissipation of flux, or equivalently, quantum irreversibility occurs for 
|A|>1. 

Following ^3], we associate to the Jacobi matrix J the spectral measure /i(x; k, A), 
defined by 

^1=/^ (26) 
J x - z 

Using the estimate (|20j). we find that the operator (J — Jo) is Hilbert-Schmidt. The 
main theorem in [Hj guarantees that the absolutely continuous component of fi(x; k, A) 
is supported in a(k, A)) = [2(1 - |A)|,2(1 + |A)|]. 

As long as |A| < 1, the support a(k, A)) is on the positive half-line. The limit 
z = ie y is straight forward, resulting in a real value for 7(fc, A). Hence for |A| < 1 
the flux in the x channel is conserved. 

When |A| > 1, the interval a(k, A)) includes the value 0. To compute we use 
(J3J), and get 

Jm 7 (M) = -^r7^(0,A)^0 (27) 

Where we write dfi(x; k, A) = g(x, X)dx. Thus, conservation of flux is violated, which is 
equivalent to the onset of irreversibility. This is the main result of the present work. It 
proves the fact that Xm / ~f(k, A) ^ 0, but it does not provide an explicit expression 
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for j(k, A). However, once it is computed numerically or approximately (see next 
paragraph) it can be used together with (j25j) to compute the vertex scattering matrix. 

<r bM {k,X)^-6 bM +t{k,X) (28) 

An approximate computation of the resolvent for energies in the vicinity of the elastic 
threshold (E > |) can be obtained by replacing the Jacobi matrix J(k, A) by Jo(A). It 
is has constant values along its diagonals, and an exact solution of the inhomogeneous 
equations for A > 1 read 



a n = e ina ^L= ; cos a = -\ (29) 
Vk ki A 



l(k, A) ~ ~T~T~ • (30) 
k h 



Thus, 



This is an explicit expression which can be directly used in(| 

Note that the valency v of the vertex enters only though the scaled strength 
parameter A = %2A. 

When we relax the condition E < Ei the same phenomenon will occur albeit it 
involves now the inelastic scattering channels, whose number is finite and determined 
by the total energy E. It gives no new insight, and therefore it will not be discussed 
further. 

The fact that the spectrum of the Jacobi matrix is continuous played a crucial role 
in the previous arguments. It appears in a way which is reminiscent of the discussion 
of quantum irreversibility presented in the introduction. 

Coming back to the observation that the ingoing current exceeds the total outgoing 
current, one naturally asks where can one find the rest of the probability density. This, 
and other relevant questions will be clarified in the next section. 

4. A simple model 

The star graph with v — 2 is equivalent to the Schrodinger equation 
d 2 1 / d 2 



dx 2 2 y dq 



+ q ) + A q 5{x) 



V(x,q) = E V(x,q) , (31) 
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with (x, q) G M 2 . The boundary condition at the vertex is replaced by a term which 
can be interpreted as a potential, linear in q and well localized in x. The Schrodinger 
equation describes a particle in the (x, q) plane, moving under the action of the potential 
shown in figure (|T|) . where the 5 function is replaced by a narrow Gaussian. In this 
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Figure 1. Grey level map of the potential in 

figure, the hight of the potential is indicated by the grey level, where black (white) 
marks the domains where the potential is most negative (positive). Consider now a 
particle moving in two dimensions and acted upon by this potential. Away from the 
domain of strong coupling, the potential is independent of x and quadratic in q, thus 
forming a valley along the x axis. Asymptotically, the particle moves with a constant 
speed in the x direction, while oscillating in the q direction. At the vicinity of x = 0, the 
potential changes abruptly. Here a positive ridge (assuming A > 0) along the positive 
q axis makes the q > domain inaccessible to the particle, while a down sloping valley 
pointing towards the negative q axis, may attract the particle. As long as we use a 
finite Gaussian to simulate the 5 potential, the quadratic potential in the q direction 
will take over, and will form a barrier from which the particle will be reflected back and 
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eventually the particle will find itself moving along the x axis in either the positive or 
the negative directions. The limit where the Gaussian becomes a 8 function, is delicate, 
and our computations show that it depends on the size of A. Once |A| > 1, the effective 
potential valley extends to infinity, and it is sufficient to to draw flux in the q direction 
without ever reflecting it back. Going back to the original way of looking at the problem, 
we can say that, with a finite probability, the "system" degree of freedom remains near 
the vertex at x = while the harmonic oscillator is excited in an outgoing stretching 
mode due to the strong interaction. 

Another point of view can be obtained by studying the stationary wave function 
(fT5|) at x = 0. We consider its projection on the space of functions spanned 
by the oscillator states with n > 1, which correspond to evanescent modes in the x 
direction, 

oo 

*fa°)= E^r*n(g). ( 32 ) 

n=l VKn 

We use the approximate expression (|29|) for a n , and replace the harmonic oscillator wave 

functions by their WKB approximation, 

12 cos S(n) 
Xn{q) ~ V^(2n + l-g2)i/4 ' 

S( n ) =( n + -)( arcsin , g + , Q (1 ? )* - - ] . (33) 

K J V 2 1 V y/2^+1 V^Tl 2n + l J 2) y } 

Substituting in (j22J), and using the saddle point approximation to perform the sum, we 
obtain: 

p—i^2~ tan a 

V(x = 0,q)^C- , (34) 
y2m cos a 

Where C is independent of q. This function is not square normalizable. 

The simple computation above provides some heuristic evidence in support of the 
suggestion that the interacting system possesses a new continuum channel, where the 
particle can be trapped indefinitely. This issue is discussed in great detail and rigor 
in Solomyak's paper ^2]- I n the sequel we shall present another variant of the simple 
model, where the onset of a continuum at |A| > 1 will be elucidated from yet another 
perspective. 
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Consider (J3Tj) subject to the boundary conditions 

W(x = -L 1 ,q)=V(x = L 2 ,q)=0. (35) 

This Schrodinger operator describes a particle which is confined to the interval [— L±, L 2 ] 
and a harmonic oscillator, which are coupled by an interaction which is very well 
localized at x — 0, and is linear in the oscillator amplitude. The spectrum of the 
uncoupled (A = 0) problem is pure point: 

E n>m = (n+ l -)+ ™ 2 i n = 0,l..- ; m = 1,2, ••• (36) 

We shall turn now to the case A ^ 0. We consider in particular the domain 
of energies E < |, just under the lowest energy of the unperturbed system. Again 
we expand fy(x, q) in the complete ortho normal basis Xn{q), and use the notation 



k n = +\/ \ + n — E for n = 0, 1, 



2 

oo 



= ~7f= , (37) 

„=0 V fc n 

where r) n (x) satisfy the boundary conditions (jSSJ, are continuous at the origin, and 



normalized by t] n (0) = 1. Explicitly, 



sinh fc n (Li+a:) £ < 
sinhfenii — 



sinh k„(L 2 -x) x > 



sinhfc n L2 

These functions are evanescent away from the origin x — 0. Substituting (|37jl in (|31jl. 
we get a set of equations from which the coefficients to be determined. 



7na„ + A ( a n+H l— + a n _ 1 /_ =0, (38) 

where, A = ^= and 

7 n = coth fc n Li + coth /c n L 2 . (39) 

Equation (|38|) is a second order recursion relation with the boundary condition 
a_i = 0. The coefficients depend on the energy E, and the spectrum of the Schrodinger 
operator is determined as the values of E for which the series a n (E)/y/k^ is square 
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summable. Before we proceed to analyze this further, we should point out that in the 
limit n » E, 



o(-) ; 7 n = l + o(-) • 



Thus, (I38|) limits to the recursion relation with constant coefficients (|21|). This is of 
crucial importance to our discussion. 

We shall now provide the heuristic arguments which show that for |A| < 1 the 
spectrum of the Schrodinger operator (}3Tf retains its pure point character. However, 
for |A| > 1 it has a purely continuous component. Using the initial conditions 
a_i = 0, ao = 1, one can apply the exact recurrence relation (|38|) and obtain cin(E) for 
arbitrary large N. Chose N such that for n > N one introduces only a small error if 
one replaces (J58J) by (J2U- The solution of (J2Tj) can be written explicitly: 

a a n pprox = A C + B C n ; f = ^ (-1 + VT^A*) , (41) 

where A, B are arbitrary constants, to be determined by the initial conditions. 

For |A| < 1, |^| < 1. Thus, to get a converging solution B must vanish. The 
matching of the approximate series at n = N requires 

a N (E) 

which is the spectral secular equation. Limiting our attention to the spectral interval 
E < |, and with L 1; L 2 sufficiently large, we can chose N = 0, and the secular equation 
reads, 

{ = ? 4a E)) )\ (43) 



oq(E) A 



Thus, 



E 



+ < i. (44) 



The approximate solution predicts the existence of a bound state of the system below 
the ground state of the unperturbed harmonic oscillator. This is rigorously proved in 
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For |A| > 1, £ = exp(— % arctan(v / A 2 — 1) = e %c > . Now, both terms in (|41j) 
contribute, and the secular equation reads, 

a N+1 (E) _ A e-« N+1 K + B eg+^ 

a N (E) ~ i e -^ + 5e^ ' ( j 

In contrast with the previous case, for every E < 1, the matching condition are satisfied 
when 

± = _ e HN + i K M^L . (46 ) 

a N (E) 

Hence, the secular equation can be solved for every E, and the spectrum is continuous. 
Moreover, the a n (E) do not vanish in the limit of large n. Substituting the a n in 
(137)1 . the resulting ty(x,q) is bounded but not square integrable, which is typical for 
eigenfunctions in the continuous spectrum. We emphasize again the fact that the onset 
of irreversibility coincides with the appearance of a continuum. This result, which was 
illustrated above for the simple model, is general, and is entirely due to the typical form 
of the boundary conditions at the vertices. 

A last comment is in order. Often, when facing an infinite set of equations of the 
type studied here, it is tempting to truncate them at a physically or computationally 
motivated point. In the present context, one would naturally attempt to truncate the 
space of closed channels, and perform the computations within the subspace of opened 
channel, exclusively. This would always result in a unitary description of the vertex 
scattering matrix, and the possibility of flux dissipation would be overlooked. 

To summarize - the main lesson from this work is that genuine quantum flux 
dissipation and irreversibility can be induced by a single harmonic oscillator - the 
transmitted and the reflected fluxes do not add up to the incoming flux. This result 
is of general interest because it shows that in a strongly interacting system, continuum 
channels, which are absent in the uncoupled system and bath, can be opened by the 
interaction, and induce dissipation. At such instances, the perturbative approach to 
dissipation might not reflect the true nature of the problem. 
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